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Abstract: It is widely believed that existing electroweak data requires a Standard Model 
Higgs to be light while electroweak and flavour physics constraints require other scalars 
charged under the Standard Model gauge couplings to be heavy. We analyze the robust- 
ness of these beliefs within a general scalar sector and find both to be incorrect, provided 
that the scalar sector approximately preserves custodial symmetry and minimal flavour viola- 
tion (MFV). We demonstrate this by considering the phenomenology of the Standard Model 
supplemented by a scalar having SU C (3) xSU L (2) xU Y (l) quantum numbers (8, 2)iy 2 — which 
has been argued [13] to be the only kind of exotic flavour singlet scalar allowed by MFV that 
couples to quarks. We examine constraints coming from electroweak precision data, direct 
production from LEPII and the Tevatron, and from flavour physics, and find that the obser- 
vations allow both the Standard Model Higgs and the new scalars to be simultaneously light 
- with masses ~ 100 GeV, and in some cases lighter. The discovery of such light coloured 
scalars could be a compelling possibility for early LHC runs, due to their large production 
cross section, a ~ 100 pb. But the observations equally allow all the scalars to be heavy 
(including the Higgs), with masses ~ ITeV, with the presence of the new scalars removing 
the light-Higgs preference that normally emerges from fits to the electroweak precision data. 
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1. Introduction 

Most physicists believe that new physics beyond the Standard Model (SM) awaits discovery 
at the LHC, and experiments at the Large Hadron Collider (LHC) will soon probe the weak 
scale and (hopefully) reveal the nature of whatever new physics lies beyond the Standard 
Model. Since the Higgs sector is among the least understood in the SM, new scalar physics 
could well be what is found. 
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However, to be found at the Tevatron or the LHC, any such new scalar physics should be 
associated with a comparatively low scale, A ~ TeV. And because the scale is low, it must 
be checked that the new physics cannot contribute to processes that are well-measured and 
agree well with the SM, such as electroweak precision data (EWPD) and flavour-changing 
neutral currents (FCNCs). This suggests taking most seriously those kinds of new physics 
that suppress such contributions in a natural way. This can be elegantly accomplished if 
the effective field theory (EFT) appropriate to low energies obeys approximate symmetries, 
such as a custodial symmetry SU(2)c [1, 2, 3] for EWPD and the principle of minimal flavor 
violation (MFV) [4, 5, 6, 7, 8, 9], which suppresses FCNCs when formulated appropriately 
[10, 11, 12]. 

Recently, it was discovered [13] that there are comparatively few kinds of exotic scalars 
that are flavour singlets and can have Yukawa couplings with SM fermions in a way that is 
consistent with MFV. The only two possible scalar representations allowed are those of the SM 
Higgs or octet scalars, respectively transforming under the gauge group SU(3) x SU(2) x U(l) 
as (l,2) 1/2 or (8,2) 1/2 . 

In this paper we examine what constraints EWPD 1 , flavour physics, and direct production 
constraints place on this general scalar sector consistent with MFV. To this end we consider 
the Manohar-Wise model, for which only one (1, 2)]/ 2 scalar and one (8, 2)^ / 2 scalar are 
present. 

Since it is the quality of SM fits to electroweak precision data that at present provide our 
only direct evidence for the existence of the SM Higgs, it is perhaps not surprising that the 
existence of a scalar octet can alter the Higgs properties to which such fits point. In particular, 
the best-fit value of the Higgs mass obtained from SM fits to EWPD is now 96±|| GeV [14]. 
We find that for the Manohar-Wise model, EWPD fits both change the implications for the 
Higgs mass, and limit the allowed mass range of the extended scalar sector. 

We find that when the masses of the Higgs and octet scalars are approximately degenerate, 
the electroweak fits allow both the Higgs and the octet to be light, with masses ~ 100 GeV 
(or even lighter for some components). Alternatively, agreement with EWPD also allows the 
octet and the Higgs doublets to be both heavy, with masses ~ 1 TeV. The Higgs doublet can 
be heavy and remain consistent with precision fits because its contribution to the relevant 
observables is partially cancelled by the contribution of the octet doublet. Having such a 
heavy Higgs without ruining electroweak fits is attractive, as a resolution of the so-called 
'LEP Paradox' [15]. We find that the precision electroweak fits generically prefer to limit the 
splittings among some of the octet components, but by an amount that does not require fine 
tuning of parameters in the potential. (The overall masses of the two multiplets are subject 
to the usual issues associated with the electroweak hierarchy.) 

The plan of this paper is as follows, in Section 2 we review the Manohar-Wise model, and 
describe its motivation as a general scalar sector that can both allow an approximate custodial 
symmetry and satisfy MFV. In Section 3 we present our results for the phenomenology of 

1 We thank J. Erler for private communication on the recent update to the EWPD fit results related to [14]. 



- 2 - 



the model. In particular, we describe its implications for an EWPD fit, and explore the 
parameter space that allows both doublets to be either light or heavy. Since the fits prefer 
a scalar spectrum that is approximately custodially symmetric, we also study loop-induced 
SU(2)c breaking, and demonstrate that the allowed parameter space is not fine tuned. This 
section also describes direct-production constraints on the Higgs and octet scalar, coming 
from both LEP2 and the Tevatron, and reexamines how previously studied flavour constraints 
change if the new octets are comparatively light. We find that the octets can pass all these 
tests, for parameters with scalars that are either light or heavy. Some conclusions are briefly 
summarized in Section 4. 

2. Theory 

In this section we recap the main features of the the model, obtained by supplementing 
the SM with an colour-octet, 5C/ i (2)-doublet scalar. Particular attention is spent on its 
approximate symmetries, since these underly the motivation to naturally satisfy FCNC and 
EWPD constraints. 

Motivation for (8,2)2/2 scalars. 

Minimal Flavour Violation (MFV) is a framework for having flavour-dependent masses with- 
out introducing unwanted flavour changing neutral currents (FCNCs). It assumes all breaking 
of the underlying approximate SU(3);y x SU(3)£> x SU(3)q flavour symmetry of the SM is 
proportional to the up- or down-quark Yukawa matrices. The fact that only scalars trans- 
forming as (8,2)2/2, or as the SM Higgs [13], can Yukawa couple to SM fermions consistent 
with MFV is the motivation of the phenomenological study we present here. 

However, we also note that octet scalars appear in many specific new- physics scenarios, 
including various SUSY constructions [16, 17], topcolour models [18], and models with extra 
dimensions [19, 20]. Various approaches to grand unification also have light colour octet 
scalars, including Pati-Salam unification [21] and SU(5) unification [22, 23, 24]. Colour octet 
doublets have also recently been used to study new mechanisms for neutrino mass generation 
[25]. Octet scalar doublets appear naturally in models of the Chiral-Colour [26, 27] type where 
QCD originates in the chiral colour group SUl(3) x SUr(3), since in this case octet doublets 
are expected in addition to the Higgs as 3 ® 3 = 8 © 1. As discussed in [28] one can also 
consider the class of models where the SM is extended with SU(N) x SU(3)c x SU(2)l x U(1)y 
and imagine model-building composite Higgs models with a (8, 2)1/2 scalar in the low energy 
spectrum. We emphasize that although many BSM scenarios contain (8,2)2/2 scalars our 
motivation is essentially phenomenological. 

2.1 The Manohar-Wise model 

In the Manohar Wise model [13], the scalar sector of the SM is supplemented with the (8, 2) 1 / 2 
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scalar denoted 



S A+ 
S A ° 



(2.1) 



where A is the colour index. 

The Yukawa couplings of the (8, 2)j / 2 scalar to quarks is determined up to overall complex 
constants, r)\j and t]d, to be 

L = nu g% u%T A (S A f e Q{ - rj D gg tf R T A {S A ^ Q{ + h.c, (2.2) 

where g u and g D are the standard model Yukawa matrices, i,j are flavor indices and 

(2.3) 



1 
-1 



The most general renormalizable potential [13] is 

,2\ 2 



V = - (H^Hi - V -\ + 2m| TV + Ai#t*iZj Tr (sft'Sj) + AaF 1 ^ Tr (s^S, 

X 3 H^H^ j Tr (SiS^ + XiH^Tr (s^SjSi) + A 5 # t4 Tr (s^SiSj) + h.c. 



+ A 6 Tr [St'SiSft'SA +A 7 Tr (.S ; .S' ; .S' '.S',j + A 8 Tr (S^J Tr (sfi'Sj 
+ AgTr (^VS,-) Tr (,S tj S 4 ) + Ai Tr (S^) Tr (V^) + AnTr (.S .V, 5 ' .V ' ) , (2.4) 



where i and j are SU(2) indices and S = S A T A . Since a field redefinition can be used to 
make A3 real, this represents 14 real parameters in the potential beyond those of the SM, 
which reduce to 9 in the custodial SU(2) symmetric case — see eqs. (2.9) through (2.12), 
below. No new parameters enter in the couplings of the (8,2)]/ 2 scalar to the electroweak 
gauge bosons since it has the same electroweak quantum numbers as the Higgs. We use this 
fact to bound the masses of the octets in Section 3.1.1. The Ai ; 2,3 terms in Eq.(2.4) lift the 
mass degeneracy of the octet states when the Higgs acquires a vacuum expectation value. 
Expanding the neutral scalar octet as 



?A0 _ D R + lD l 



S AV = ' J (2.5) 
V2 



the tree level masses become [13] 



v 2 

Af| = Mj + Ai — 



Mr = Ml + (Ai + A 2 + 2A 3 ) — 



Mf = Mf + (Ai + A 2 - 2A 3 ) j. (2.6) 
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2.1.1 Custodial symmetry 

We find below that EWPD fits prefer the masses of some of the scalars in these models 
to be approximately degenerate in mass. In particular, fits prefer a mass pattern that can 
be naturally understood as being due to an approximate custodial SU(2)c symmetry, under 
which the SM vector bosons transform as a triplet and the Higgs transforms as a singlet and 
a triplet. This symmetry is broken in the SM both by hypercharge gauge interactions, and 
by the mass splittings within fermion electroweak doublets. 

For these reasons we next explore the implications of the custodial-invariant limit, for 
which SU(2)c is an exact symmetry of the underlying new physics beyond the SM. In this 
scenario, it is interesting to examine the case that SU(2)c is preserved in the Manohar-Wise 
model potential at a high scale ~ 1 TeV, up to the breaking that must be induced by the 
SM. Imposing exact SU(2)c on the octet Higgs potential we find that the potential can be 
rewritten in terms of bi-doublets 



*=(e^*,$, S A = (eS\,S A ), 



(2.7) 



where e is given in Eqn. (2.3) and the most general gauge- and custodial-invariant potential 
becomes 



V 



Tr [$t$ 



16 

+ai Tr (cS f $) Tr (V $) + (bi Tr[T A T B T c 



+ -Tr Tr 



$ f $ ) Tr ( S A J S A 



]Tr [&S A S] J Sc)+h.c 



+ciTr[T A T B T c ]Tr [S\Sc) Tr(5^$ 



+cZ 1 Tr[T A T B T c T D ]Tr 
+ei Tr[T A T B ] Tr[T c T D ]Tr 
+/i Tr[T A T B ] Tr[T c T D ]Tr 



S\S B 



S a Sb) TtIS^Sd 



S\S C ) Tr(5' S D 



:t 



(2.8) 



where Ta is used as a basis in colour space with 9 independent terms when the potential 
is SU(2)c invariant. 2 Expanding out the potential and comparing to the general result of 
eq. (2.4), we confirm the result of [13] that SU(2)c implies 



2A 3 -- 
2A 6 : 
A 9 = 

but we also find the additional constraint" 



A 2 , (2.9) 
2A 7 = An, (2-10) 
Aio, (2.11) 

A 4 = A*. (2.12) 

2 An alternative way to obtain this count is to regard SU(2)l xSU(2)c as 50(4), with both H and S A 
transforming as real fields in the 4-dimensional representation. In this case the invariants of the potential can 
be written m 2 s (S A ■ S A ), d ABC {H ■ S A ){S B ■ S c ), f AB a{Hi ■ SfS% -Sf) f kl ,(H-H)(S A ■ S A ), (H ■ S A )(H ■ S A ), 
(S A ■ S A ) 2 and the two independent ways of colour-contracting (S A ■ S B )(S C ■ S D ). 

3 We thank A Manohar for communication on this point clearing up a subtlety. 
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Note that this constraint can effect the production mechanism of the octets at Tevatron and 
LHC. We see in particular that because SU(2)c symmetry implies A2 = 2A3, in this limit M± 
and Mi become degenerate. 

2.2 Naturalness issues 

In general, even if the scalar potential is required to be custodial invariant at a particular 
scale, it does not remain so under renormalization due to the presence of custodial-breaking 
interactions within the SM itself. In this section we compute these one-loop symmetry break- 
ing effects, allowing us to quantify the extent to which the custodial-invariant potential is 
fine-tuned. To do so we calculate in Feynman gauge and note that ghost fields do not couple 
to the components of the S doublet. We also neglect goldstone boson contributions to the 
mass splitting as they come from the SU(2)c symmetric potential and so therefore cancel out 
in the mass splittings; not leading to mixing between the Sr and Si states. 

SU(2)c breaking due to Yukawa corrections 

The breaking of SU(2)c due to Yukawa couplings is straightforward, the requisite diagrams 
are given by Fig 1. 



s * \ s- s> v & s< s < 

b a) t b) "~' b " V 

I b 

& \^ J s H s' V J s R 

Figure 1: SU(2) violating contributions to S 1 , S masses from the yukawa sector of the theory. 

The correction to the mass S~ S + two point function comes from diagram (a) and is 
given by 

5(T[S + S-}) Y = 3 b — — + — !^S£M^j + A °( m t } ~ P 2B o(P 2 > m l m t)\ / 213 x 

16 TT 2 V 2 

(ml \i] D \ 2 + m\ \r) V \ 2 + ml m 2 (\r] D \ 2 + |r?[/| 2 - 2 ijp rju - 2 rf D r^)) B {p 2 ,m 2 b , m 2 )) 

lb 7r z f; z 

where we express our results in terms of Passarino-Veltman (PV) functions whose definitions 
are given in [42], and we set |Vtb| — 1- 

The contributions to the S] operator comes from the diagrams (b) and (c) and is given 

by 

x/tiqIqI\\ x rn 2 t (2A (m 2 )\i lu \ 2 + B (p 2 , m 2 ,m 2 )(A m 2 Im [rju] 2 - p 2 \ r]u \ 2 ) ) 
d(T{S S \) Y = -dab ' 
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ml(2A (ml)\r] D \ 2 + B (p 2 ,ml,ml)(4mllm[r) D } 2 -p 2 \r] D \ 2 )) 
- o ab .(2.14) 

We are interested in the mass splitting of Mf and M± , however to the accuracy we work 
one can also easily calculate the shifts to 5{T{S R S R })y and 5(T{S R S I })y due to the mixing 
induced between the real and imaginary components of S A0 . With these results we can then 
obtain the contributions to the diagonalized Mj. The correction to 5(T{S R S R })y is given 
by the same diagrams as 5(T{S I S I })y with the appropriate replacements, giving 

R R m 2 t {2A {m 2 t ) \r] V \ 2 + B (p 2 , m 2 , m 2 t ){4m 2 t Re[rm} 2 - p 2 \r]u\ 2 )) 
6{T{S S })y = -d ab 16^2 ' 

m 2 (2A (ml) \r] D \ 2 + B (p 2 , m 2 , m 2 )(4 m 2 Re[ry D ] 2 - p 2 \yp\ 2 )) 

16 ir z v z 

The mixing of the Sr, Sj fields at one loop 5(T{S R S I })y is given by diagrams (d,e) and 
is given by 

x/rrfaR a i\\ x {mlRe[r] D }Im[ri D } B (p 2 ,m 2 b ,m 2 b ) - m\ Re[rj u ] Im[w] A)0 2 , m t, ™t)) 

6{T{S S }}y = -dab ^2^2 

which is only nonzero when at least one of the MFV proportionality constants rjD,r]u are 
imaginary as expected. We define the mixing angle and renormalize the theory in the Ap- 
pendix. 

Gauge sector SU(2)c violating corrections 

Calculating the required four diagrams represented by diagrams (g,i) in Fig 2 one finds 



W^» G = ^_8 AB (^MMM + dMMM _ h A p^M 2 w ^M 2 ± \ - ^-W,MlMl 



where cw = cos[#iy] and the integral is given in terms of PV functions as follows 

I 3 [p 2 , M 2 ,M 2 ] = {2p 2 + 2M 2 - M 2 )B [p 2 , M 2 , M 2 ] + 2A [M 2 ] - A [M 2 ]. (2.16) 

The result for 5(T{S R S r })g is identical up to the replacement Mr — > Mi. One can 
similarly calculate the other six diagrams corresponding to (f,h) that give the following con- 
tribution for 5(T{S + S~})g in terms of PV functions 4 

,2 \2 



s (r{s + s - })o _ _g_ { - (?MMM + iStzlMMM& _ \ hb \Mi,Ml\ 

-\h\^',Ml„M'j] - (1 " 4 ^ >2 J3lp 2 ,Ml, Ml) - s'h\p\ 0, Mi]) (2.17) 



,2 ^,2 
3. 

u 

Mixing between the states S 1 , S R is forbidden in the gauge sector as the couplings are real. 



4 Note that diagram (f) with a photon loop is scaleless and vanishes in dim reg. 
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rJVV. w.z,a rS^h w - z 



_ L<aaI zhyyS 



s + s s> s> 

h) i) 

Figure 2: SU(2) violating contributions from the gauge sector of the theory. 



Given these loop-generated effects, we wish to estimate how large the custodial-symmetry- 
breaking interactions are once we run down to observable energies from the scale of new 
physics. The answer depends on how far we must run, however due to the hierarchy problem 
of the Higgs mass (which is only accentuated when more light scalars are added to the 
spectrum), it is likely that new physics must intervene at a relatively low scale for new 
physics of ~ TeV. Such a low scale for a UV completion implies that the symmetry structure 
of the UV theory is consistent with EWPD and flavour constraints. 

The splitting induced by SM interactions is given by the difference between the renor- 
malized mass at A and the low scale, where we ignore the running for simplicity in this 
estimate 

f' m 8 M 2 

/ (^)d^ = Mi[Z^ = A)-Z^ = m)}, (2.18) 

where Zf di is the leading perturbative correction of the mass counterterms, whose values are 
given explicitly in the appendix using a zero-momentum subtraction scheme. 

As is shown in detail in the next section, the largest Mj,M± SU(2)c violating mass- 
splitting that is allowed by our EWPD fit is approximately ~ 40(55) GeV for the entire 
68%(95%) confidence regions (see Figure 6). We now examine how natural such a small 
splitting is assuming a typical low mass of 150 GeV. 

In determining the splitting, the values of rji employed are critical. For the lower bound 
on the r\i we take the one approximate loop radiatively induced value rji ~ 0.35 2 /(167r 2 ). 
Note that we use the result of [40] that determined an upper bound on \r]u\ from the effect of 
the octet on = [Z — > bb)/(Z — > Hadrons). For charged scalar masses of (75, 100,200) GeV 
the one sigma allowed upper value for \r]u\ is (0.27,0.28,0.33). 

For M± = 150 GeV, we choose the couplings to give the largest induced splitting con- 
sistent with other experimental constraints (rju = 0.3, rjo = 0.45), Mj = 150 GeV (its value 
before the perturbative correction in the high scale SU(2)c preserving scenario) and Mr = 
(190, 230) GeV which are the maximum values consistent with EWPD for the (68%, 95%) re- 
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gions. We find that the EWPD regions begin to have tuning for a high scale degenerate mass 
spectrum at (90 TeV, 8000 TeV). Conversely choosing the unknown rju,rjD ~ 0.35 2 /(16 ir 2 ) 
one finds that the (68%, 95%) regions begin to have some degree of tuning for scales of 
(170 TeV, 19000 TeV). For a UV completion that approximately preserves MFV and SU(2) C , 
considering a SM and octet low energy scalar mass spectrum allowed by EWPD is not a fine 
tuned scenario. 

3. Phenomenology 

We next turn to the various observational constraints. As we shall see, the most robust 
constraints are those coming from the absence of direct pair-production at LEP, which require 

M± £ 100 GeV and M R + Mj£ 200 GeV . (3.1) 

Since the octet scalar couples to both photons and gluons, these constraints are essentially 
kinematic up to the highest energies probed by LEP (more about which below). 

3.1 Fits to Electroweak Precision Data 

A strong restriction on the properties of exotic scalars comes from precision electroweak 
measurements, whose implications we now explore in some detail. The dominant way that 
such scalars influence the electroweak observables is through their contributions to the gauge 
boson vacuum polarizations; the so-called 'oblique' corrections [31, 32, 33]. The calculation 
of the oblique corrections proceeds as usual with the vacuum polarizations being determined 
directly by evaluating the diagrams given in Figure 3. 

^\ 

vAAAi k/WV vAAAAAA/^ 

a \ y b a b 

Figure 3: Self energies calculated for the EWPD constraints on the octets, where (a, b) = 
(W + W~ , ZZ, 77, Z'y). The self energies needed to determine STUVWX are given in the Appendix. 

When evaluating these it is important to keep in mind that the direct production con- 
straints, eq. (3.1), can allow one of Mr or Mj to be significantly lower than 100 GeV. This is 
important because it precludes our using the most commonly-used three-parameter (S, T and 
U) parametrization of the oblique corrections [31, 32, 33], since these are based on expanding 
the gauge boson vacuum energies out to quadratic order: U a b(q 2 ) ~ A a b + B a f,q 2 , where a 
and b denote one of Z, W or 7. Since the electroweak precision measurements take place at 
q 2 ~ or q 2 ~ M|, using the quadratic approximation for U a b(q 2 ) amounts to neglecting 
contributions that are of relative order M|/M 2 , where M is the scale associated with the 
new physics of interest (in our case the new-scalar masses). This approximation becomes in- 
adequate for M below 100 GeV, and so we must instead use the full 6-parameter description 
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(STUVWX), such as in the formalism of ref. [29, 30]. In general, the STUVWX formalism 
reduces to the three-parameter STU case when all new particles become very heavy. 

For ease of comparison with past results we start by quoting the results we obtain for the 
fit to the six parameters of the STUVWX oblique formalism, regardless of how they depend on 
the parameters of the Manohar-Wise model. The results are given in Table 1, which compares 
the results obtained by fitting 34 observables (listed in an appendix) to (i) all six parameters 
(STUVWX); (ii) only three parameters (STU); or just two parameters (ST). The number of 
degrees of freedom in these fits to (6,3,2) parameters is v = (28,31,32), respectively. The 
X 2 /v for the three fits is within one standard deviation \Jljv = (0.27, 0.25, 0.25) of the mean 
of 1, indicating a good quality of fit. The experimental values and theoretical predictions 
used are given in Table 2 in the Appendix. 



Oblique 


STUVWX Fit ( X 2 /v = 0.91) 


STU Fit ( X 2 /v = 0.99) 


ST Fit ( X 2 /v = 0.98) 


S 


0.07 ±0.41 


-0.02 ±0.08 


-9.9 x 10" 3 ±0.08 


T 


-0.40 ±0.28 


-0.02 ±0.08 


1.1 x 10~ 2 ±0.07 


U 


0.65 ±0.33 


0.06 ±0.10 




V 


0.43 ±0.29 






w 


3.0 ±2.5 






X 


-0.17 ±0.15 







Table 1: EWPD Fit Results in various schemes for the 34 observables listed in the Appendix. The 
STU and ST fits fix the other oblique corrections to zero as a prior input. The error listed is the 
square root of the diagonal element of the determined covariance matrix. The central values of the 
fitted oblique corrections decrease as more parameters are turned off. All three fits arc consistent with 
past results and the PDG quoted fit results. 



The correlation coefficient matrix for the three fit results are as follows, 



Mstuvwx 



1 

0.60 
0.38 
-0.57 


-0.86 



0.60 
1 

-0.49 



0.38 
-0.49 
1 



-0.95 0.46 
-0.01 
-0.13 -0.76 



-0.57 
-0.95 
0.46 

1 


0.13 



-0.86 \ 
-0.13 
-0.01 -0.76 

0.13 

1 

1 / 



(3.2) 



M S TU 



1 0.84 
0.84 1 
-0.20 -0.49 




M 



ST 



1 0.87 
0.87 1 



(3.3) 



We use the results of this fit to constrain the masses allowed in the Manohar-Wise model 
by computing the vacuum polarizations as functions of the masses of the octet and Higgs 
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3i> r ' f '- 

50 100 150 :>' 250 300 

,W+ 



Figure 4: Comparison of the three and six parameter fits for low masses. (The upper two panels 
are not symmetric about Mj = Mr and Mr = M + because we scan only through positive values for 
the couplings, A;.) The three parameter fit is red (grey) and the six parameter fit is blue (black). 
Contrary to naive expectations the six parameter fit is more constraining on the model despite the 
extra parameters; the correlations between the extra parameters (S,X and U, X and T,V) increases 
the constraints on the model. The masses are in GeV. EWPD constrains the mass spectrum to be 
approximately SU(2)c symmetric in cither case where M± ~ Mj. 



scalars. We obtain allowed mass ranges for the scalars by demanding that the contribution 
of the new physics (and the difference between the floating Higgs mass and its fiducial value, 
which we take from the SM best fits to be 96 GeV), A% 2 which satisfies 

(C -1 )ij(A0i) (A9j) < 7.0385 (12.592) (3.4) 

for the 68% (95%) confidence regions defined by the cumulative distribution function for the 
six parameter fit. Here C is the covariance matrix constructed from the correlation coefficient 
matrix given in eq. (3.2) or (3.3) 

x _i aV(fl) 



-n- 



and A Oi = Ai - A{ u is the difference in Ai = S, T, U, V, W, X as a function of octet masses 
and the best fit value, given in Table 1. 



Y 




AX once marginalized over Y 



Figure 5: A cartoon of the best-fit confidence interval for a strongly correlated pair of variables, 
indicating how the best constraints can be missed once one of the variables is marginalized. 

An example of the best-fit regions for the allowed octet masses is given in Figure 4, 
which compares the quality of the constraints that are obtained using the full six-parameter 
(STUVWX) parametrization, as opposed to the three-parameter (STU) expression. The three 
panels plot the masses of the components of the octet that lie within the 68% confidence 
ellipsoid of the best-fit value as the various scalar couplings, Aj, are varied. The two panels 
of this plot show how these masses are correlated by the condition that the predictions agree 
with the precision electroweak measurements, and the points in the upper two panels all 
satisfy Mj < Mr and M + < Mr because we choose to scan only through positive values of 
the couplings Aj. 

The strongest correlation is between Mi and M+, for which agreement with EWPD 
demands these two masses cannot be split by more than about 50 GeV. This is as might be 
expected given that this difference must vanish in the limit that the potential is custodial 
invariant. The breaking of SU(2)c generically leads to bad fits because custodial-breaking 
quantities like the parameter p — 1 = a T are measured to be very small. New sources of 
SU(2)c breaking are extremely constrained, the deviation of the p parameter from its SM 
value is given by p = 1.0004+°;™°* [39]. 

The comparison in Figure 4 also shows that the six-parameter STUVWX fit agrees with 
the three-parameter STU fit when all scalars are heavy, as might be expected. It also shows 
that the six-parameter fit is the more constraining one when the octet masses are light. 
We understand that this happens because of the strong correlations amongst the oblique 
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parameters, which implies that the best-constrained parameter direction is not aligned along 
any of the STUVWX clXCS , clS shown in Figure 5. As a result the constraint obtained by 
restricting to the axes V = W = X = can be weaker than the full result, significantly 
affecting the determined 68% confidence regions. For this reason our remaining results quote 
only the results of the full six-parameter fit. 




50 100 150 200 250 



Figure 6: Comparison of the 68% red (grey) and 95% blue (black) confidence regions when < \ < 1. 
The masses are in GcV, and Mj,M + < Mr because we scan only through positive values of the 
couplings A,; . For low masses the 95% confidence region is significantly expanded compared to the 68% 
region, this is due to the spread of available masses being larger for low masses, as the mass splitting 
between the states scales as ~ v 2 /m s . We examine the naturalness of this mass spectrum in Section 
2.2 and find that it is not simply a fine tuned solution for an underlying new physics sector. 

3.1.1 Constraints on Octet scalars 

Figure 6 displays the 68% and 95% confidence regions of the model for couplings that range 
through the values < A» < 1, while Figure 7 does the same for couplings that run through 
the larger range < Aj < 10, where i = 1,2,3. As noted above, agreement with the 
EWPD selects an approximately SU(2)c symmetric mass spectrum, where A2 ~ 2 A3 and 
\M± — M/| < 50GeV, but this is easily understood. Consider the case where the octets are 
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heavy, v 2 /Mg -C 1, which was examined in [13]. In this mass regime it is the model that 
constrains the mass spectrum to be degenerate, M± ~ Mr ~ Mj, since the mass splittings 
scale as v 2 /Ms from Eq. (2.6). The contribution of the octets to the S and T parameters, 5 is 
then [13] 

A v 2 

s= Mii> T = 96^^p^; ( ^- (2A3)2) ' (3 - 6) 

where sw = sin(#jy)- Large corrections to S and T are avoided if Aj decreases and preserves 
approximate SU(2)c as Mg decreases, therefore allowing smaller octet masses. 




Figure 7: Comparison of the 68% red (grey) and 95% blue (black) confidence regions when A, < 10. 
Notice that the region selected for by EWPD for Mj ss M + that is approximately SU(2)c symmetric 
is not enlarged. 

How natural are the small intra-octet splittings favoured by EWPD? If the mass splitting 
is induced by the potential, while v 3> M s , for the octet masses to be allowed by EWPD that 
selects for a mass degeneracy A M = Mj — M± , one would have to require that the couplings 
the the octet-Higgs potential satisfy the scaling rule 

A 2 - 2A 3 < 4 TAT. (3.7) 

v 

5 We have checked that our results in the STUVWX formalism reduce to these results when u 2 /M| <C 1. 
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As EWPD requires A M ~ 50 GeV for the 95 % confidence region this is a mild hierarchy of 
couplings given by A2 — 2 A3 <C 0.8 \f\\. Conversely for the case mg 3> v, one requires that 
the couplings the the octet-Higgs potential satisfy the scaling rule 



which is easily satisfied for small Aj (which we see below are favoured by Landau pole con- 
straints). 

The calculations presented in previous sections for the running of custodial-breaking 
couplings can be used to frame a criteria as to whether the above coupling pattern is natural. 
The scale dependence of the masses is used to estimate what the SU(2)c splitting of the 
masses should be in the theory below the UV scale, A, without tuning. One determines 
how high the scale A can be before the EWPD mass regions are excluded. This quantifies 
the degree of fine tuning of the masses for this scenario. 6 Since the electroweak hierarchy 
problem argues that the scale of new physics is likely not too much larger than the TeV 
regime, we find that the favoured mass splittings are natural, provided that the underlying 
theory approximately preserves MFV and SU(2)o 

The above ranges of allowed splittings amongst scalar masses directly constrain the three 
couplings Ai 2,3 to be small. But small Aj, for i > 4, are also favoured due to considerations 
of the effect of these Aj on the running of the Higgs self coupling [28]. The mild assumption 
that one not encounter a Landau pole while running the Higgs self coupling up to 10 TeV, 
when one assumes Aj>4 = and m/, = 120 GeV, gives the constraints [28] 



However, generically Aj>4 7^ and if the octets and the Higgs were part of a new sector 
then the cut-off scale could be lower that 10 TeV. For these reasons we only take these 
constraints to inspire the Aj < 1 limit for the parameter space searches in Figure 6, but also 
examine parameter space where we relax this bound to Aj < 10 in Figure 7. We emphasize 
that direct production bounds on the octets that rely on their fermionic decays essentially 
constrain the MFV proportionality factors rji, while EWPD is complementary in that it 
constrains the parameters in the potential, Aj, by constraining the mass spectrum. 

3.1.2 Implications for the inferred Higgs mass 

Adding the new octet scalar to the SM also affects the best-fit value of the Higgs mass that 
emerges from fits to EWPD. In particular, we now show that the presence of the octet can 
remove the preference of the data for a light Higgs, even if the new octet scalar is also heavy. 

6 To determine the mass splitting, we technically need to diagonalize the Si field which mixes at one loop 
with S R . As the non diagonal terms in the mass matrix are one loop, the effects of this diagonalization on 
the mass eigenstate Sj shifts the mass at two loop order. See the Appendix for a determination of the mixing 
angle. Thus to one loop order one can just take the one loop corrections to Mi and M± of the last two sections, 
properly renormalized, to determine the mass splitting through the counterterms. 



A 2 - 2A 3 < 8 



(AM) m s 

v' 2 





(3.9) 
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To determine this effect we calculate the one-loop Higgs contribution to the six oblique 
parameters and jointly constrain the Higgs mass and the octet masses in the fit. For example, 
S in this case becomes 

S = S oct {M R , Mj, M±) + S mggs {M h ) - S Higgs (M h = 96 GeV) (3.10) 

where S oct ( H iggs) is the one-loop octet (Higgs) contribution to the S parameter. We neglect 
the two-loop dependence on the Higgs mass in the fit and this leads to an underestimate of 
the allowed parameter space, as we find the 68% (95%) confidence level values of fitting the 
Higgs mass alone are given by 112 (160)GeV. This gives a conservative range when comparing 
to the various allowed values that are strongly dependent on the priors used in the PDG . 

The effect of the octets changes the preferred Higgs mass significantly, and two mecha- 
nisms are at work depending on the size of the octet mass. If the octet mass Mj is small, it 
can allow the Higgs mass to increase by effectively replacing it in the oblique loops, thereby 
giving agreement with EWPD. This is illustrated in the upper-left plot of Figure 8, which 
shows how a large Higgs mass correlates with small Mj. 

The other panels of Figure 8 reveal another mechanism at work, however 7 . In these one 
sees that as the upper limit on Aj is increased, the upper limit on the Higgs mass confidence 
regions becomes significantly relaxed. This is due to a cancellation between the effects of the 
heavy octet and the Higgs in their contributions to oblique parameters, that is made possible 
by a positive AT contribution that the octets give to x 2 - For the three-parameter fit, the x 2 
test is of the form 

(C-%-(A0i) (A9j) = 596 (AS) 2 - 1159 {AS) {AT) + 751 (AT) 2 (3.11) 

where we neglect contributions that are not logarithmically sensitive to the Higgs mass at 
one loop, since this is all that is relevant to the argument. For the three-parameter fit, the 
68% confidence region is defined by (C~ )jj(A0j) {A6j) < 3.536 and is easily satisfied for 
light Higgs masses. As the Higgs mass grows, its contribution to {AS) and (AT) becomes 
dominated by the logarithmic dependence 

< A *>^(3) - d (AT, -iM3)' (3 - i2) 

where Mh is the reference value of the Higgs mass, which for our fit is 96 GeV. The crucial 
point is that (AT) is negative for Mh > Mh and for the SM this quickly excludes large Higgs 
masses because of the sign flip in the {AS) (AT) term in \ 2 . 

Including the contribution of the octets in the large mass regime (v 2 /Mj <C 1) modifies 
these expressions to 



12vr \ M H J 6vrM 2 ' 



7 Note that we expect a careful study of the non oblique Higgs and octet mass dependence of Rb will further 
constrain this parameter space with all scalars heavy but not remove it. 
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Figure 8: The effect of octets on the fitted value of the Higgs mass. The plots of Mh versus the 
other octet states are substantially the same. The green line is the 68% confidence bound where the 
Higgs alone is varied at one loop. The yellow line is the 95% confidence bound where the Higgs alone 
is varied at one loop, and the black line is the direct production bound on the Higgs mass at 95% 
confidence. The red (grey) region is the 68% confidence region, while the blue (black) region is the 
(95%) confidence region for a joint fit to the octets and the Higgs. Notice the increase in vertical scale 
for the diagrams as the upper limit of the Ai is increased through 1 (upper left), 3 (upper right), 6 
(lower left) and 10 (lower right). The mechanism that is allowing the Higgs mass to increase and still 
be in agreement with EWPD is the postitive AT contribution from the octets that is discussed in 
Section 3.1.2. 



<^> - -m* (§) + 9e^W Ai - (2A3)2) ' <3 - 13) 

where the factor X 2 — (2 A3) 2 comes from a factor of (M| - M|)(M| - M|) in the octet 
contribution, and is a measure of the total mass splitting in the doublet. For Aj > 0, we 
know Mp > M± and so the octets give a positive contribution to (AT) so long as M? > M±. 
The octets (or any other doublet with gauge couplings and small mass splittings) then allow 
(AT) in Eqn. 3.13 to be positive, and so allow a large degree of cancellation between the 
(AS) 2 , (AT) 2 and (AS 1 ) (AT) terms in Eqn. 3.11. The size of the positive (AT) contribution 
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scales with the upper limit on Aj , explaining the significant relaxation of the Higgs mass bound 
in Figure 8. We find that the Higgs and the octet scalars could both have masses ~ 1 TeV 
and still lie within the 95% contour mass region allowed by EWPD. We also note that we 
restrict our searches to positive Aj (which must be so for at least some of the couplings to 
ensure the absence of runaway directions in the potential) , however clearly negative A2 could 
also act to relax the EWPD bound on the Higgs mass by giving a negative contribution to 
(AS). 

We emphasize the generic nature of the mechanism, wherein the contributions of TeV 
scale new physics can mask the contributions of a heavy Higgs to electroweak precision ob- 
servables. It applies in particular when EW symmetry breaking leads to a mass splitting 
of an extra SU(2) doublet, since the extra doublet can give a positive contribution to (AT) 
proportional to the mass splittings of the doublet components. This has been recognized as 
a simple way to raise the EWPD bound on the Higgs mass by satisfying the positive (AT) 
criteria of [35]. Expressed as an effect on the p parameter, it also has a long history going 
back to observations by Veltman [34], being rediscovered for two-Higgs-doublet models in 
[36], and used for the construction of the Inert Two Higgs doublet (IDM) model [37]. 8 In this 
latter model, the Higgs mass is raised, addressing the "LEP paradox", and the naturalness 
of the SM Higgs sector is also improved by raising the cutoff scale of the modified SM. In the 
IDM model a parity symmetry is imposed to avoid FCNC's. 

We note that the example of the general scalar sector consistent with flavour constraints, 
the Manohar-Wise model examined in this paper, naturally has a number of the benefits 
of models like the IDM while avoiding the imposition of a parity symmetry. Allowing the 
second doublet to couple to quarks improves its potential for detection, without introducing 
large FCNCs due to MFV. It is interesting that the effect of raising the Higgs mass has 
emerged naturally from the most general MFV scalar sector and was not a model building 
motivation of the MW model. Variants of the MW model, can address the naturalness of 
the scalar sector through raising the cut off scale and further the colour charge of the octet 
provided some rational for the second doublet not obtaining a vev, through the avoidance 
of the spontaneous breaking of colour. Also, for the entire parameter range, octets skew the 
distribution of the allowed Higgs masses so that the direct production bound on the Higgs 
mass and the EWPD fit of the Higgs mass can be in better agreement. 

3.1.3 Implications for the tension between leptonic and hadronic asymmetries 

Although the SM produces a good quality global fit to EWPD, there exists a mild tension in 
the data between the leptonic and hadronic asymmetries. In particular A b FB deviates from the 
SM predicition by 2.5a and favours a heavy Higgs ~ 400GeV, while A e differs from the SM 
by ~ 2a and favours a Higgs mass far below the direct production bound. Here we address 
the question of whether the oblique contributions of octet scalars can change this tension. 

To this end we calculate \ 2 f° r the hadronic asymmetries A b FB , A C FB , A^, A c , and for the 
leptonic asymmetries using A T and the A e values given in Table A. The results are shown 
8 For a similar construction see [38] 
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Figure 9: The \ 2 °f the leptonic asymmetries (red) and hadronic asymmetries (blue) as a function 
of Higgs mass in GeV. The solid curves show the contribution of the Higgs alone and the dashed 
curves are for the Higgs and the octets. The figure on the left is for octet masses (M±, Mr, Mj) = 
(300,400,330)GeV and on the right is for (M±,M R ,M I ) = (900, 1000, 940) GeV. 



in Figure 9, where the solid curves plot \ 2 with the SM Higgs alone and the dashed curves 
include the octets for a particular mass spectrum allowed by EWPD. The two panels compare 
results for relatively light and relatively heavy octet scalars. 

The figure shows that the preferred value of the Higgs mass is strongly dependent on 
the mass splitting of the octets. As discussed in Section 3.1.2, the octets, unlike the Higgs, 
give a positive contribution to AT, which depends on the mass splitting in the doublet. This 
increases the allowed value of the Higgs mass. The octets can change the pull of A e , for 
example, to favour large Higgs masses, however they also do the same to A b FB . As can be 
seen from Figure 9, although the leptonic and hadronic asymmetries can now both prefer a 
Higgs masses above the direct production bound of 114.4 GeV, they are not brought in to 
closer agreement in their predictions for the value of Mh- 

We see from this that the octet oblique contributions do not in themselves remove the 
tension between the leptonic and hadronic asymmetries. However, because the octets are 
coloured it is possible that their non-oblique corrections to A FB might be able to bring 
together the leptonic and hadonic observables. We leave this observation to a more complete 
calculation, which lies beyond the scope of this paper. 

3.2 Direct-production constraints from LEP 

The octets would have been directly produced at LEP2 if they were light enough through the 
processes in Figure 10. 

The production cross sections are given by 




A 3 / 2 (l,— ±,— ± ) (3.14) 



MJ | iM z Vz 



i 2 / 2 , 2\ 

+ v + (v e + a e ) 
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Figure 10: The tree level production mechanism for S + + S and + Sj at LEPII. 



d A /4vra 2 \ 3/2 / Ml M 2 \ , 2 , Af| iM z r z 2 



(3.15) 



where we have defined = 8, a e = — (4swc\y) 1 

X(x : y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz, (3.16) 

s w~ c w 1 — /„ ,m 
f + = — — , vq = , v e = — — (3-17) 

2s\yC\Y 2s\yC\v 4:S\yCW 

The highest COM energy at which LEP2 operated was y/s = 209 GeV, where approximately 

J Cdt ~ 0.1 fb _1 of integrated luminosity was collected. We give a rough estimate of the 

sensitivity of LEP2 to light octets by requiring less than 10 total events for a given set of 

masses, a x J Cdt < 10. Note that these limits are essentially kinematic limits for production, 

and more accurate exclusions in the mass parameter space are possible, but these will be 
dependent on the detailed decays of the octets and SM backgrounds and be weaker constraints. 
The LEP2 production bounds are shown in Figure 11. 

3.3 Tevatron constraints 

3.3.1 Dijet constraints on the production cross section. 

Heavy octet production via gluon fusion has been examined in some detail in the literature see 
[13, 40, 27]. We use the results of [13, 40, 27] to determine the production cross sections for 
light octets and consider the relevant bounds on the model in this region from the Tevatron. 
The single production cross section we use, [40], neglects for simplicity the scalar mass splitting 
and assumes that r]u, A4 and A5 are real. However, note that this is partially justified for light 
masses as EWPD selects for an approximately degenerate mass spectrum with an approximate 
SU(2)c symmetry in the underlying potential, giving A4 = A£ and one need only assume one 
of the couplings are real. 9 For the sake of simplicity we will also neglect the effects of mixing 
of the Si, Sr states that can occur if the effective yukawa couplings of the octet carries a 
phase as discussed in the Appendix. The pair production cross section for the charged scalars 
is twice that for the real scalars [13] and so is not shown. 



9 Note that setting A4 and A5 to real values removes the scalar loop contributions to the single production 
of Si, which can become large as the values of A^s increases. 
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The tree level pair production dominates the loop suppressed single production in the 
low mass region for small A45. However as A^s increase the single production contribution 
takes over, which occurs at A^s ~ 2 for the neutral scalar, Sr, with a mass of 200 GeV. 

A direct search strategy to find octets is to look for narrow resonance structures above 
the QCD background for states that decay into dijets. CDF has recently performed such a 
search [45] with 1.13 fb~ 1 of data that could discover octet bound states [41] or single Si that 
decay to dijets above the QCD background. The cross sections for the production of these 
states at the Tevatron, leading to dijet resonance structures, are orders of magnitude below 
the QCD background in the regions of parameter space we consider, this is shown in Fig. 12 

The low mass region is not directly ruled out, although a dedicated study to refine the 
lower mass bound is warranted due to the shape dependence of the exclusion bound. 10 



Other possible indirect search strategies for the effects of octet scalars include determining the effect of 
the octets on the Af b . In a similar manner to axigluons [47], these new exotic coloured states could contribute 
to A'ft, as they are coloured, couple strongly to tops, and are not a vectorlike state. Interestingly, Af b has 
recently been measured [48, 49] to be A FB = 0.19 ± 0.065(stat) ± 0.024(syst) which is a deviation larger than 
2 sigma from its SM value [47] of A FB = 0.05 ± 0.015. 
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Figure 11: Comparison of the 68% (red or light) and 95% (blue or dark) confidence regions when 
Xi < 1. The LEP2 production bound for ten events is the black line. 
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Figure 12: Shown is the production cross section of a(g g — > red short dashed line, c(gg — + 5/) 
blue long dashed line, and the cr(gg — > SrSr) given by the solid green line. The results are for 
Tevatron with ,/i = 1.96 TeV, a s (M z ) = 0.1217, m t = 173.1 GeV, M z = 91.1876 GeV and the NLO 
CTEQ5 pdfs. The values of (A4, A5) chosen are (0,0) upper left, (1,1) upper right and (10,10) for 
the bottom graph. In all three graphs we have set T]u = 0.2. The dependence on r]xj is weak and as 
T)u decreases the production cross sections decrease. Also shown is a 95% confidence limit band (the 
shaded region) derived from [45] that places an upper bound on new physics that decays to dijcts. The 
region is defined by the upper limit on <r(X)B(X — > jj) * A(\y\ < 1) where the difference between the 
W and RS graviton G* 95% confidence upper bounds are taken and the acceptance fraction requires 
the leading jets to have rapidity magnitude \y\ < 1. The exclusion region depends weakly on the 
shape of the resonance, so a dedicated study is required to exactly bound the octet decay to dijets, 
however, the octet signal is orders of magnitude below the exclusion regions obtained from Tevatron 
before branching and acceptance ratios further reduce the signal. A rcsummation of large threshold 
logarithms for single S production was performed in [27]. The K factors for single S production was 
found to be ~ 2 for 500 GeV a octet mass and this K factor falls as the mass decreases. This indicates 
that threshold enhancements will not raise the cross section enough to exclude octets for the entire 
low mass region. 

3.3.2 Gauge boson decays and Lepton Signatures 

The decays of the octets involving gauge bosons 

S RJ ^W ± ST, Srj^ZS^r 
S ± -> W ± S RJ , S ± -» ZS ± . (3.18) 

were studied in some detail in [13, 28]. These decays are of phenomenological interest as the 
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gauge bosons can be a source of leptons to trigger on at LHC and Tevatron. The EWPD 
constraints \M± — Mj\ < 50GeV and for most of the allowed parameter space \Mi — MA < 
Mw,Mz, as the mass splitting of the doublets scale as v 2 /M s for large masses. This causes 
the decays to proceed through an offshell gauge boson for most of the allowed parameter 
space. In this case an effective local operator can be used to approximate the decays. 

For example consider Sr — > S~ t + v through an off shell W . The effective Lagrangian at 
leading order is given by the product of scalar octet and left handed lepton currents 

C eS = ^L(S R d fl S + )(u Ll »£ L ). (3.19) 
V 2 A% 

Exact formula for three body decays such as this exist in the literature [43]. For the masses 
allowed by EWPD 11 generally the energy release is A = Mr — M± < Mr,M_,Mw- The 
resulting decay width at leading order in A /Mr is 

a 2 A 5 , , 

(3.20) 



60vr^M^ 



When Mr > 2 mt the decays to leptons through an offshell W, Z are suppressed decay 
channels. The dominant decay widths are to tb, ti unless rju <C t)d- The ratio of Ti to this 
decay, in the limit Mr > 2mj, is given by 



O.OOSGeV/ A 



s o^ t t M R \rnj\ 2 V 50GeV 



for a = 1/128, s w = 0.48 and m t = 173 GeV. 

When Mr < 2 mt the offshell W, Z will be dominant decay channels for light masses for 
much of the parameter space. Taking nib = 4.23 GeV, and the other factors as before, the 
ratio of the offshell decay to the S R — > bb decay is given by 

T f 4a 2 ( A 5 v 2 



8 GeV / A \ 



M R \r] D \ 2 V 50 GeV 



5 



1 . (3.22) 



If the dominant fermionic decays are to charm quarks due to a mild hierarchy of rjjj > 
(m&/m c ) rju, then taking m c = 1.3 GeV gives the branching ratio 

82 GeV ( A 

(3.23) 



> c5 M R \ m \ 2 V50GeV 



lx This assumes that the initial state that is eventually triggered on is not highly boosted. This is generally 
the case due to the kinematic reach of the Tevatron and LHC. 



- 23 - 



Thus when quark decays are suppressed through Mr < 2rri( the dominant decay mode 
will be through an offshell W, Z for much of the parameter space of r]u,VD allowed by other 
constraints, notably the constraints due to Rj,. This sets a lower bound on the decay width 
of the heavier octet species given parametrically by Eqn. 3.20. This sets an upper bound 
on the lifetime of these components of the octet doublet of 4.5/A 5 ps which yields a upper 
bound on the decay length of the form 1CP 3 /A 5 m. 12 Thus the heavier octet species will 
decay promptly inside the detector and not leave a long lived charged track signature. 

As dominant decay modes of the heavy components of the octet doublet (when Mj < 2 mt) 
can be three body decays, the final state signature would be excess monojet or dijet (depending 
on the boost of the final state octet) events in association with a lepton and missing energy, 
or enhancements of dilepton signatures with a monojet or dijet. Dedicated studies of these 
signatures are warranted. The lifetime of the lightest component of the octet doublet is 
dictated by its decay to fermion pairs. 

3.3.3 Constraints from tt decays. 

For neutral octet masses above 2 mt, decays into top quark pairs can be dominant. These were 
previously considered in [40]. The observed limits on excess ax • B{X — > tF) at Tevatron with 
0.9 fb~ l of data [46] do not rule out octets in the intermediate mass region 350 — 1000 GeV. 
The production cross section for single gg — > Sr production can become large enough for 
the bound on tt to be relevant, however this requires A4 ~ A5 ~ 75 which is well into a 
nonperturbative region of the potential making any conclusion suspect. We illustrate these 
limits in Fig. 13 

3.3.4 Constraints from bbbb decays. 

The dominant decays for light masses will be to quarks 5+ — > tb, Srj — > bb below the ti 
threshold for r]u,D ~ ®(X)- m this regime [28] places a lower bound on the scalar mass of 
approximately 200 GeV from the CDF search for a scalar particle decaying dominantly to 
bb when produced in association with b quarks [44] This bound is avoided for almost all of 
the available parameter space for light octet masses. S^r can decay preferentially to charms, 
which corresponds to a mild hierarchy of couplings 

1 770 1 2 rn 2 . 1 

rnk < -I ~ t?^ 3 - 24 
\vu\ m b 10 

when neglecting 0{m 2 c /Mg) terms. Neutral scalar masses below 200GeV are allowed for 
Vd ^5 0.1, given an upper limit of rju ~ 0.3 from [40] for masses in this range. The three 
body decays discussed in Section 3.3.2 are actually dominant over quark decays for much of 
the parameter space allowed by EWPD for light octet masses, invalidating the assumptions 
of [28] for most of the remaining parameter space. 



12 Here we have converted units assuming that A is given in GeV as a pure number, ie for A = 50 GeV we 
have a upper bound on the lifetime of 1.2 x 10~ 2 as. 
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Figure 13: Shown is the production cross section of a(g g — > Sr) red short dashed line, a(g g — > 5/) 
blue long dashed line, and the <r(gg — ► SrSr) given by the solid green line. The results are for 
Tevatron with y/s = 1.96 TeV, a s (Af z ) = 0.1217, m t = 173.1 GeV, M z = 91.1876 GeV and the NLO 
CTEQ5 pdfs are used. The DO 95% confidence limit on a(X)T(X — ► tT) is the upper solid black 
line [46]. The values of (A^As) are (10,10) for the left hand figure and (75,75) for the right hand 
figure. r]u = 0.2 for both figures. For perturbative A, < 10, current Tevatron production bounds on 
resonances in it do not rule out octets of mass 350 — 1000 GeV. 



3.3.5 Constraints from 77 decays. 

A promising signature for octets at hadron colliders is the annihilation of a pair of charged 
octets to photons, gg — ► S + S~ — ► 77. We can use the recent results of DO [51, 50] that 
utilize 4.2 fb^ 1 of date to place 95% confidence upper limits on cr(h) x BR(h — ► 77) compared 
to the SM Higgs signal to directly constrain octet annihilation into 77. We must consider 
annihilation decays of octet bound states, octetonia, studied in [41], as the contribution from 
virtual octets will be a non-resonant signal and the Tevatron Higgs search would not apply. 
Due to the fact that the results are reported only up to Higgs masses of 150 GeV we are 
only able to exclude octets up to 75 GeV, which is already disfavoured by LEP2. If the 
experimental study of h — > 77 is extended to higher Higgs masses at the Tevatron or LHC, 
this signal is likely to be a significant constraint on the model. 

We utilize the fact that this signature has been studied for octetonia in [41] to demonstrate 
the potential of this signal to raise the mass limit on octets. The ratio we are interested in is 
that of the octetonia a(gg — > + ) x BR(0 + — > 77) to the SM rate for a(gg — > h) X BR(h — > 
77). We take [41] 

a(gg - + ) x BR(0 + - 7 7) - ^^^'V " m 2 /s) (3.25) 

where s is the partonic center of mass energy squared and |^(0)| is the wavefunction at the 
origin. We have used the approximation BR(0 + — ► 77) ~ a 2 /«g(2Ms). For the Higgs, we 
take the approximation 
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a{gg -+ h) x BR{h - 77) « ^ ( ) 10 " 3 ^ " M h/*) ( 3 - 26 ) 
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Neglecting order one factors the ratio of these two signals scales as 

« Ml *£Mp!( *<L) (3 , 7) 

This ratio must be less than ~ 35 [51, 50] for M h = (100, 150) GeV or M± = (50, 75) GeV. 
Unless the wavefunction at the origin was much smaller than its approximate expected value 
given by [41] 

. , . , l9 N^a 3 Jmsv)M^ 

h/>(0)| 2 = ' ( 3 - 28 ) 

O 7T 

this bound will likely be violated for this entire mass range. Extending this analysis to higher 
Higgs masses is expected to raise the lower mass bound on octet states for this reason. For a 
recent comprehensive study of octetonia signals in gamma gamma for octets from ~ 200 — 500 
GeV see [41]. 

3.4 Flavour constraints reexamined for light scalars 

Flavour constaints on (8, 2) 1 / 2 scalars were examined in some detail in linear MFV 13 in [13] 
when the masses of the octet scalars were considered to be ~ TeV. However, although MFV 
suppresses flavour changing effects and ensures the vanishing of tree level flavour changing 
neutral currents in linear MFV, when one goes beyond leading order in the Yukawa couplings 
problematic flavour changing neutral currents are possible [40]. The correct way to examine 
such flavour issues is to utilize a nonlinear representation of MFV 14 such as formulated in 
[10, 11, 12] which is beyond the scope of this work. 

We have reexamined the flavour constraints that were examined in [13] in linear MFV 
for the light octet masses allowed by EWPD and not ruled out by direct production bounds. 
Flavour constraints are largely irrelevant for \rjjj\ once the far more restrictive constraint from 
Ri) is known. To quantitatively demonstrate this consider K° — K° mixing for relatively light 
masses M s = 300 (400)GeV. We use the results of [13] for the contribution of the octets to the 
wilson coefficient (C s ) of the operator (V* d Vt s )(dL "1 V slX^l j u s£) and use the SM expression 
of [52] for the contribution of this operator to K° — K° mixing and hence \tx\- One finds 
that the contribution of the octets to | ex | is given by 

A\e K \ = \C e B K lm[V? d V ts ] Re[Vf d V ts ] C s \ (3.29) 

Using the measured values m K = 497.6 MeV, f K = (156.1 ± 0.8)MeV, (AM K ) exp = 3.483 ± 
0.006) x 10 _12 MeV one obtains 

C e = F * — ^ = 3.65 x 10 4 . (3.30 
6V2tt 2 AM k 



13 Where one only utilizes a linear yukawa coupling for the scalars. 
14 We thanks J. Zupan for discussions on this point. 
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Further, Lattice QCD [53] gives the input i?^(2GeV) = 0.54 ± 0.05, and using the central 
values of fitted values for the CKM parameters A, fj, p, A from the PDG we find that the shift 
in | | is given by 

A|ejr| = 1.5 (1.6) x 10- 12 (| W | 2 + 6 (3)| W | 4 ) (3.31) 

for M s = 300 (400)GeV. Considering \e K \ exp = (2.229 ± 0.010) x 10~ 3 while the same values 
employed above gives the central value \ex\theory = 1-70 x 10~ 3 one can set an upper limit on 
\r)u\ from K° — K° mixing by conservatively assigning one tenth of the discrepency between 
theory and experiment to the effect of octets. This gives an upper bound on \rju\ of 48 (56) 
for M s = 300 (400)GeV. The weak mass dependence of the bound allows one to neglect 
Kaon mixing constraints for low masses, compared to Rb constraints on \r]jj\, for light masses 
M s < 1 TeV, in linear MFV. 

The B — > X s 7 decay rate constrains the combination ItjuWdI, m the limit rjjj is small, 
and was calculated in [13] . Using their result and the upper bound on \rju\ from Rj,, we 
determine the strongest upper bound on \t]d\ f° r light masses by requiring that the octet 
contribution to B — > X s -y is less than the ~ 10% SM theoretical and experimental errors. 
For M± = (75, 100, 200) and the corresponding maximum \rju\ = (0.26, 0.27, 0.33), one obtains 
an upper bound on \rjjj\ of (0.36,0.39,0.50). As \rju\ decreases, the upper bound on \rj£>\ is 
relaxed. 

Finally, the electric dipole moment of the neutron constrains the imaginary part of the 
iji and using [13] we find for light masses that Im[r?^ rf^\ < 1/10 for mj = 100 GeV. 

4. Conclusions 

We have considered the phenomenological constraints of the general scalar sector that contains 
one (1,2)^2 Higgs doublet and a one (8, 2)u 2 colour octet scalar doublet. To this end we 
have performed a modern fit in the STU and STUVWX approaches to EWPD and used 
these results to determine the allowed masses for light octets. We have demonstrated that, 
somewhat surprisingly, the six parameter fit formalism is more restrictive for light states due 
to strong correlations amongst the fit observables. We find that the octet doublet masses can 
be in the 100 GeV range. Such light octets can significantly effect the discovery strategies for 
a light Higgs by modifying the Higgs production mechanism through a one loop contribution 
to gg — > h that is not well approximated by a local operator. Octets will also induce a further 
effective coupling at one loop between h and 77, ZZ and W + W~ and would significantly 
effect Higgs discovery at LHC [54]. Despite this, we have shown that current production 
bounds on light octets at LEP2 and Tevatron do not rule out the low mass region and further 
studies for narrow resonances in the dijet invariant mass distribution and h — > 77 signal 
are required. Currently, octets are another example of physics beyond the SM that can 
significantly effect the properties of the Higgs and yet are otherwise relatively unconstrained 
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experimentally. For light octets, one possible alternate search strategy is to utilize the Higgs 
Pt distribution [57] to find indirect evidence for onshell octet scalars that have eluded direct 
detection. 

We have also performed a joint fit for the Higgs and the octets by varying the Higgs 
mass oblique corrections at one loop while allowing the masses of the octets to vary. Doing 
so we have demonstrated a mechanism that is quite general in its effect of giving a positive 
contribution to the T parameter when an extra doublet is present and fit to in EWPD. This 
allows the Higgs and octet to be simultaneously heavy and the Higgs can be as massive as 
its unitarity bound. For the parameter space where the Higgs mass is raised, h decaying to 
pairs of octets is kinematically suppressed. The search strategy for the heavy Higgs remains 
substantially the same with difficulties in constructing a mass peak due to the width of the 
Higgs resonance and large irreducible backgrounds to SM processes producing W + W~ decays 
such as from it, and large Wj backgrounds. Likewise very heavy octets are also very broad 
resonances for large masses and are difficult to discover at hadron colliders with decays to it 
dominating, and large SM backgrounds. Further dedicated studies of LHC phenomenology 
of this scenario are warranted, as are further dedicated studies to attempt to raise the lower 
mass bounds on octet scalar doublets. 
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A. EWPD fit 

The data and theory predictions used in constructing the fit are given in Table 2. 



Observable 


Data Used 


Theory Prediction 


M w [GeV] 


80.428 ± 0.039 


80.380 ± 0.015 




80.376 ± 0.033 


80.380 ± 0.015 


M z [GeV] 


91.1876 ± 0.0021 


91.1874 ± 0.0021 


T z [GeV] 


2.4952 ± 0.0023 


2.4954 ± 0.0009 


T had [GeV] 


1.7444 ± 0.0020 


1.7419 ± 0.0009 


T inv [MeV 


499.0 ± 1.5 


501.68 ± 0.07 


T l+l - [MeV] 


83.984 ± 0.086 


84.002 ± 0.016 


<Jhad [nb] 


41.541 ± 0.037 


41.483 ± 0.008 


i? e 


20.804 ± 0.050 


20.736 ± 0.010 


Rfj, 


20.785 ± 0.033 


20.736 ± 0.010 


R T 


20.764 ± 0.045 


20.736 ± 0.010 


R b 


0.21629 ± 


0.00066 


0.21578 


± 


0.00005 


R c 


0.1721 ± 


0.0030 


0.17224 


± 


0.00003 


A e 

^FB 


0.0145 ± 


0.0025 


0.01627 


± 


0.00023 


AH 

FB 


0.0169 ± 


0.0013 


0.01627 


± 


0.00023 


A T 
FB 


0.0188 ± 


0.0017 


0.01627 


± 


0.00023 


A b 
FB 


0.0992 ± 


0.0016 


0.1033 


± 


0.0007 


A c 
FB 


0.0707 ± 


0.0035 


0.0738 


± 


0.0006 




0.2316 ± 


0.0018 


0.2315 


± 


0.0001 


A e 


0.15138 ± 


0.00216 


0.1473 


± 


0.0010 




0.1544 ± 


0.0060 


0.1473 


± 


0.0010 




0.1498 ± 


0.0049 


0.1473 


± 


0.0010 


\ 


0.142 ± 


0.015 


0.1473 


± 


0.0010 


A T 


0.136 ± 


0.015 


0.1473 


± 


0.0010 




0.1439 ± 


0.0043 


0.1473 


± 


0.0010 


A b 


0.923 ± 


0.020 


0.9347 


± 


0.0001 


A c 


0.670 ± 


0.027 


0.6679 


± 


0.0004 


9l 


0.3010 ± 


0.0015 


0.3039 


± 


0.0002 


& 


0.0308 ± 


0.0011 


0.03000 


± 


0.00003 


9 v v e 


-0.040 ± 0.015 


-0.0397 ± 0.0003 


97 


-0.507 ± 0.014 


-0.5064 ± 0.0001 


Qw(Cs) 


-73.16 ± 0.35 


-73.16 ± 0.03 


Qw(Tl) 


-116.4 ± 3.6 


-116.8 ± 0.04 


T w [GeV] 


2.141 ± 0.041 


2.0902 ± 0.0009 



Table 2: Observables used in fit to oblique parameters. 
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The numbers we use for the theory predictions are based on the 2008 PDG results of a 
global fit to the EWPD. The input values used in the theory predictions are 



M z = 91.1876 ± 0.0021GeV, M H 
m t = 173.1 ± 1.4GeV, a s {M z ) 
a{M z y l = 127.909 ± 0.0019, A cnjj^ 



96l^GeV, 

0.1217 ±0.0017GeV, 

0.02799 ±0.00014. 



(A.l) 



The definitions of the oblique corrections we use are 



aS 



4s 



aT 
aU 



4,2 



aX 



Sn zz (M 2 z ) - 5U ZZ (0) 



Ml 



2 

W 



Ml' 



sw cw 



8U' Z J0)-5U' (0) 



5U WW (0) 5U ZZ (0) 



M 2 W 



Ml 



5U WW (M W ) 



«I W (0) 



4 



c w 



6IL zz (Ml) - 5n zz (0) 



Ml 



sw OT' 77 (0) -2s w cw Sn.' Z y{0), 



aV = 6IL' zz (Ml) - 
aW = SU! WW {M W ) 



5U ZZ (M 2 ) - 5U ZZ (0) 



Ml 



5U WW (M W ) - 5U WW (0) 



-sw cw 



5U Zl (M 2 ) 
Ml 



«4 7 (0) 



(A.2) 



The self energies to determine these results are given by the following in terms of PV 
functions that match the definitions in [42] and are 



16 7rV" n 



16 7rV" n 



16 7rV" n 



16 7rV" n 



d n > 



1 



d n > 



i (2 7r) n q 2 — m 2 + i e 
1 



i(2 ir) n 


[q 2 


- mf + i e] [(q - 


p) 2 


— m 2 , + ie] 


d n q 










i(2-n) n 


[q 2 


— m\ + i e] [(q — 


p) 2 


— m?, + ie] 


d n q 










i(2 ir) n 


[q 2 


— m\ + ie) [(q - 


P) 2 


— m| + ie] 



A (m 2 ) (A.3) 

B {p 2 ,ml,ml) 

P^Bi(p 2 ,ml,ml) 

P^iPv B 2 i{p 2 ,ml,m 2 ,), 
+9t,uB 2 2(p 2 ,ml,m 2 ] ) 



Our results are 



5U ww (p 



2, 9i r 



B 22 (p 2 , Ml Ml) + B 22 {p 2 ,Mi, M\ 



■ 

2ir 2 

\mm%) - \mm%) - \a (m 2 ) 



(A.4) 
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SU zz (p 2 ) 



9i 



(1 - 2s 2 w ) 2 ( B 22 (p 2 ,Ml,Ml) - ^A (M 2 ' 



+ B 22 (p 2 ,M 2 R ,M 2 ) - \A Q {M 2 R ) - \a (M 2 ) 



5U^(p 2 ) 



2e 2 



7T^ 



B 22 {p 2 ,M 2 + ,M 2 + )- l -A {M 2 + ) 



5U, z (p 2 ) = e9l{1 - 2s ^ 



TT 2 C\V 



B 22 (p 2 ,Ml,Ml)--A (M 2 



For p 2 = these expressions become 



5U WW (0) = |L Q/(M+,M B ) + i/(M +) Mj) 



OT zz (0) 



where 



./'•;"'.• "'2 i '"f • "'2 2 2 2 



9 2mfm, , mf 
- 1 — " log ,', 



The derivatives of the vacuum polarizations are 

OT; 7 (0) = -£^B (0,M 2 + ,Ml 



e gi (l-2s 2 w , , , 



12tt 2 c\y 



-Bo(0,M±,Mi 



(A.5) 
(A.6) 
(A.7) 

(A.8) 
(A.9) 

(A.10) 

(A.ll) 
(A.12) 



B. Renormalization 



We use dim reg in d = 4 — 2 e dimensions. We introduce wavefunction renormalization and 
mass renormalization constants for the octet fields as usual 



St 



s 



(o) 



Mi 



M 



(o) 



(B.l) 



However, in choosing renormalization conditions, we note that to define the masses and the 
mass splittings one cannot use MS, as in MS the mass is defined to have only the divergence 
subtracted from the bare mass. The resulting renormalized mass in MS is not shifted by 
the finite components of the loop corrections that we have determined. The renormalization 
prescription we use is the zero-momentum subtraction scheme [58], where we require that the 
self energy and its derivative with respect to external momentum, p 2 , vanishes at p 2 — > 0. 
Note that the second derivative term in the Taylor expansion of the self energy does not 
contribute until two loop order and therefore can be neglected here. The counter terms in 
the lagrangian are given by 



]T [(Zi - s> d^Si) - (Zt z Ml - \)m 2 sf] 



(B.2) 
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With this prescription the wavefunction renormalization and the mass counterterm are of the 
form 



Z t = l 

e 

b 



a d^i(p 2 ) 



dp 2 lp ^° 

„2\ 



Zmi = 1 + - + ) | p 2^ +(1 - Zi 



(B.3) 



where a, b are the coefficients of the p 2 , M 2 dependent one loop divergences respectively and 
the Sj are the finite terms of the one loop self energy. 

Using this scheme and the divergence properties of the PV functions, the wavefunction 
renormalization factors are determined to be 



9l 



7 1 yilm/r + yblvpl , 

1 64vr 2 e 32^ 2 e 

yflm[r]u] 2 +y 2 lm[r] D } 2 g| 
48tt 2 16vr 2 



1 + 



1 



2r 2 



+ 



2/t 2 lwl 2 log ^ +v 2 \vd\ 2 log 



32 vr 2 



Mo,<,n4) + MO -i-^ 

1 c w 



9\ 



7 = 1 _ ^Iwl+^M 
fl "~ 64vr 2 e ~ ' 32^ 2 e 

y 2 Re[ W ] 2 + j/ 2 Re[7 ?1) ] 2 5l 2 



1 + 



2r 2 



+ 



2/t 2 |w| 2 log ^ + y 2 |rf log 



32 vr 2 



48 ^ 2 



16 7T 2 



6o(0,M^Ml) + -^ M - M ^ 



7 , i/t 2 l^l 2 + yglrf , s 2 [ 1 , (i-2^) 2 2 - 

Z± = 1 64^ + 32^ 1 + _ ~ 2 + 2Sw 



1e 2 



1c 2 



+ 



9 i, , 9 9n (1 - 2s 2 v ) 2 b (0,Ml,Ml) , 
6o(0, M^, M|) + 6o(0, M^, M 2 ) + ^ ^ ° 2 V ' Z ' ±J - 4s 2 w log 



32 7T 2 

0/ 2 M 2 + y 2 M 2 

32 7T 2 



-w 



b (0,m b ,m t ) 



Ml 



- 1 

(B.4) 



Using these results the mass renormalization factors are determined to be 



Zmi = (2 — Zj) 



32vr 2 M 2 ' 
+y 4 (Re[r /D ] 2 + 3Im[? ?D ] 2 



% 4 (Re[ W ] 2 + 3Im[r ?r/ ] 2 



2e 



- log 



mi 



I 1 





"m 2 " 


( log 




V2e 6 


U 2 J 



5 2 



64 vr 2 M 2 . 
32tt 2 M 2 



,,2 ,, 2 (3 Ml -Ml) 

3M^-Ml + ± z 2 r; 

— ■ c 



(M^ - 2M|) 6o[0, A%, M±] + 



*> 2 {vi\mj\ 2 + vt\riD\ 2 ) 

32vr 2 M 2 



(M 2 -2M 2 )6 [0,M z ,M i? ] 



2r 2 



+Ml I 1 - log 



MjU 1 - 2 log 



M 2 



2r 2 



1 - 2 log 



L M 2 J 



Ml 
1c 2 





\M 2 R 




(l - log 


)] 


[ V 2 _ 
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ZmR — ZmI Im|— Mj,Z I ^Z R ,Ke^lm 



Zm± = (2 - Z±) - 



2 4 I 12 
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' 64tt 2 M| 
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log 
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log 
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4M| 
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9l 



64tt 2 M| 
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The remaining renormalization is for the mixing operator Sr Sj which is renormalized 
as usual by introducing a further counter term to subtract the only divergences of composite 
operators as in MS 



Zri 



(B.7) 



where 



Zr-1 Zj-\ j/fRe^JIm^] -y 4 Re[?7 D ]Im[?7 D ] 
^ = 1 + — 2~ + -2~ + 



(B.8) 



C. Mixing of S R and S/ 

For completeness in examining one loop effects we determine the mixing between Sr and Si. 
The mass matrix is given by 



Mri 



' Mj + SiTiS 1 S I }) G + 6{T{S T 5 J }) y 5(T{S R S J })y 

5(T{5 K S J }) Y M 2 + ,5(T{5 R S*}) G + S(T{S R S R }) Y 



(C.l) 
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We diagonalize the mass matrix by introducing a mixing angle and rotating the Sr, Si fields 
to a diagonal mass basis S' R , S'j via 



Si] _ ( cos(0) sin(0) \ ( S'j 
Sr ~ I -sin(0) cos(0) ) I S' R 



(C.2) 



The mixing angle is given by 



■ to\- \ y t B^{p 2 ,mlmj)Re{ri u )lm{r lu ) - yf B^{p 2 ,ml,ml)Re(rj D )Im{ri D )\ 
sm{V) — - l^-'jj 

7T A 



2 



where .Bq is the usual PV function with the divergence subtracted given by 

B*(p 2 , ml m 2 ) = -2 + log (™p) - log (i±|) (C.4) 



where /? = yl — 4mf/p 2 , which would be the velocity of the scalar produced in the CM 
frame which was subsequently to mix into another state with mass mj. We take p 2 = m 2 
as the mass splittings are a small perturbation in a radiatively induced mixing. If we take 
fi ~ 1 TeV as the scale at which we impose exact SU(2c) on our scalar potential, this gives a 
mixing angle 

sin(^0.04 I^IN^l (C.5) 

which depends weakly on the value of m s as the numerical coefficient changes by 25% for m s 
varying between 0.01 — 300 GeV. This mixing angle, if non zero, will effect the production 
cross section of the Si, Sr states at LHC and Tevatron, and introduce mixing between the 
octetonia states discussed in [41]. 
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